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13.  ABSTRACT  (Muimum  200  mpoNEU 

The  overall  objective  of  the  research  program  is  to  develop  and  test  an 
improved  model  for  the  process  of  molecular  diffusion  in  turbulent  reactive 
flows.  In  application  to  turbulent  combustion,  a  major  shortcoming  of  ex¬ 
isting  models  is  that  they  are  non-local  in  composition.  A  model  has  been 
developed,  based  on  the  construction  of  a  Euclidean  minimum  spanning  tree 
(EMST).  This  model  is  inspired  by  the  mapping  closure,  and  reduces  to  it 
|  in  the  case  of  a  single  composition.  In  general,  the  model  is  asymptotically 

local,  and  hence  overcomes  a  major  flaw  in  previous  models. 

Additionally,  studies  have  been  made  of  stochastic  Lagrangian  models  for 
turbulent  reactive  flow's:  and  an  exact  expression  has  been  obtained  for  the 
probability  density  function  of  temperature  (or  other  random  quantities)  in 
statistically  stationary  turbulence. 
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1  INTRODUCTION 


The  treatment  of  molecular  diffusion  remains  one  of  the  major  difficulties  to  be  overcome  in 
models  of  turbulent  combustion.  One  view  of  the  overall  problem  is  obtained  by  considering 
how,  at  a  given  point  within  a  turbulent  combustion  device,  the  fluid  composition  changes 
with  time.  There  are  three  processes  that  cause  this  change:  convection,  reaction,  and  molec¬ 
ular  diffusion.  In  pdf  methods  the  first  two  of  these  processes  are  treated  exactly,  while  the 
third — molecular  diffusion — has  to  be  modelled.  Many  other  theoretical  and  experimental 
studies  lead  to  the  same  conclusion:  the  major  current  issue  in  turbulent  combustion  is  to 
understand  and  model  the  effects  of  molecular  diffusion. 

A  most  promising  recent  advance  is  the  development  of  mapping  closures  (Chen  et 
al.  1989,  Kraichnan  1990,  Pope  1991,  Gao  1991).  This  is  a  new  formalism  that  yields 
“constant-free”  pdf  closures.  In  its  initial  application  to  the  marginal  pdf  of  a  scalar  (Pope 
1991)  the  accuracy  of  the  mapping  closure  has  been  remarkable. 

For  simple  test  cases,  analytic  solutions  to  the  mapping-closure  equations  can  be  obtained. 
But  for  application  to  inhomogeneous  flows  of  practical  importance,  the  closure  needs  to  be 
implemented  as  a  particle  method.  For  a  single  scalar,  a  particle-implementation  has  been 
developed  by  Pope  (1991).  The  first  topic  reported  (in  Section  2.1)  is  an  extension  of  these 
ideas  to  the  case  of  multiple  scalars. 


In  the  first  18  months  of  the  project,  a  significant  fraction  of  the  effort  was  devoted  to 
developing  a  mapping  closure  for  the  joint  pdf  of  composition  and  its  gradient.  The  addition 
of  the  composition  gradient—  compared  to  pdf  methods  currently  in  use  provides  valuable 
information  on  the  microstructure  of  the  composition  field  that  is  responsible  for  molecular 
diffusion.  One  virtue  of  this  approach  is  that  it  can  reconcile  flameiet  and  nonflamelet 
notions  in  both  premixed  and  diffusion  flames  (Pope  1990).  For  the  joint  pdf  of  composition 
and  its  gradient  contains  the  same  information  as  the  joint  pdf  of  mixture  fraction  and  its 
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dissipation  rate,  that  is  the  cornerstone  of  flamelet  models  of  diffusion  flames  (Peters  1984). 
Similarly,  this  joint  pdf  contains  full  information  on  the  surface-to- volume  ratio,  which  is 
used  in  flamelet  models  for  premixed  combustion. 

The  approach  adopted  required  the  development  of  a  computational  procedure  to  gen¬ 
erate  a  volume  preserving  mapping  between  a  Gaussian  field  and  a  field  with  specified  joint 
pdf  of  amplitude  and  gradient.  Some  of  the  progress  made  is  described  in  last  year's  re¬ 
port.  However,  this  problem  resisted  our  concerted  efforts,  and  we  decided  to  abandon  the 
approach  in  favor  of  the  more  tractable  one  described  in  Section  2.1. 

The  process  of  molecular  mixing  is  not,  of  course,  independent  of  the  velocity  field  that 
is  convecting  and  distorting  the  composition  fields.  In  pdf  methods,  there  are  separate 
stochastic  models  for  the  evolution  of  velocity  and  composition  following  fluid  particles. 
Two  papers  concerning  these  models  and  their  interconnection  have  been  written  and  are 
outlined  In  Section  2.2. 

Finally,  a  new  exact  n  suit  has  been  obtained  for  the  pdf  of  any  stationary  random  process, 
such  as  the  temperature  at  a  point  in  a  turbulent  combustion  device.  This  is  described  in 
Section  2.3. 

The  outlook  for  the  final  year  of  the  grant  is  presented  in  Section  3. 

2  ACCOMPLISHMENTS 

2.1  EMST  Mixing  Model 

The  term  '  mixing  model”  refers  to  a  turbulence  sub-model  that  describes  the  evolution  of 
the  pdf  of  composition.  In  the  Lagrangian-pdf  framework,  the  mixing  model  specifics  how 
the  composition  </>(£)  evolves  following  a  fluid  particle.  In  one  popular  model  (IEM)  </>(£) 
relaxes  to  the  local  mean  value  (<j>)  at  a  specified  rate.  In  another  class  of  models  (particle- 
interaction  models),  the  composition  of  the  n-th  particle  in  an  ensemble  d>{n)(t)  changes  by 
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an  exchange  with  another  randomly  selected  particle  (tri,  say,  with  composition  )). 

Such  mixing  models  have  been  extensively  examined  for  inert  Hows  (e.g.  Pope  15)82) 
and  several  shortcomings  have  been  identified  and  are  now  well-appreciated.  More  rt'cently 
a  different  shortcoming — peculiar  to  reacting  flows — has  been  identified.  Specifically,  the 
physics  of  the  problem  shows  that  mixing  is  local  in  composition  space  whereas  the  models 
cited  above  are  non-local.  These  ideas  are  now  explained  in  more  detail. 

In  the  simplest  case  of  Fickian  diffusion  with  uniform  diffusivitv  T.  the  rate  of  change  of 
composition  4>  due  to  diffusion  is  simply  F V2d>.  Since  the  Laplacian  is  a  local  operator,  one 
^ouu  think  of  changing  as  a  result  of  composition  differences  in  an  infinitesimal  neighborhood 
of  the  point  in  physical  space  (x).  And  since  the  composition  field  is  mathematically  smooth, 
the  compositions  </>  at  points  in  this  neighborhood  differ  infinitesimally.  Thus  mixing  is  local 
in  both  physical  and  composition  spaces.  In  both  of  the  models  mentioned  above,  the 
composition  $n)(t)  is  influenced  non-locally,  either  by  (<£)  or  by  d>{rn)(f). 

The  EMST  mixing  model,  now  described,  is  an  asymptotically  local  model,  inspired  by 
the  mapping  closure.  For  a  single  composition  <j>  it  reduces  exactly  to  the  amplitude  mapping 
closure. 

Consider  first  a  single  scalar  0(x.t)  in  homogeneous  turbulence.  In  a  particle  method, 
the  pdf  of  d>  is  represented  by  an  ensemble  of  N  particles,  the  n-th  having  the  scalar  value 

For  this  problem,  the  mapping  closure  yields  a  fascinating  particle  method  (Pope  1991). 
Let  the  particles  be  ordered  so  that 

0w(t)  <d>m(t)  <  ...  <0{N)(t).  (1) 

Then  (according  to  the  mapping  closure)  the  particles  evolve  bv  a  coupled  set  of  ode’s 

jt4>(n)d)  =  /?n+. (d>(n+,)  -  <?(n)]  +  [0{n-])  -  0{n)],  (2) 
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where  (for  large  N)  the  positive  coefficients  are 


and  A  is  a  known  function.  (The  coefficients  B \  and  /3V,  :  are  zero,  so  oll))  and  olA  +  1' 
can  remain  undefined.)  A  simple  interpretation  of  Eq.  (2)  is  that  o[n)(t)  is  drawn  to  its  two 
neighbors  at  a  rate  proportional  to  its  separation  from  them.  Note  that  as  .V  tends  to  infinity, 
the  difference  between  and  o^n±1)  tends  to  zero.  Hence  the  method  is  (asymptotically) 
local. 

Consider  now  the  general  case  in  which  there  are  a  compositions:  in  a  typical  combustion 
problem  (employing  simplified  kinetics)  a  may  be  5.  Then  the  composition  of  the  n-th 
particle  is  denoted  by 

<*‘"’(0  =  {©!"'(<),  oSr’fo.  ■■ .  (-it 

and  it  can  be  regarded  as  a  point  in  the  (7-dimensional  composition  space.  Thus  the  ensemble 
of  N  particles  correspond  to  N  points  in  a- space. 

The  concept  of  ordering  -as  used  in  Eq.  (1) — is  peculiar  to  one-space.  There  is  no  direct 
equivalent  in  multi-dimensional  spaces,  and  so  a  direct  extension  of  the  model  defined  by 
Eq.  (2)  is  not  possible. 

However,  we  have  developed  a  model,  inspired  by  Eq.  (2),  which  is  asymptotically  local, 
and  which  reduces  to  Eq.  (2)  in  the  one-composition  case.  It  is  based  on  Euclidean  minimum 
spanning  trees  (EMST).  An  example  of  an  EMST  (for  a  =  2.  .V  =  100)  is  shown  in  Fig. 
1.  By  definition,  the  EMST  is  the  set  of  edges  joining  the  points,  such  that  all  points  (ire 
connected,  with  the  N  -  1  edges  chosen  (out  of  the  N 2  possibilities)  so  that  their  total  length 
is  minimal.  By  this  construction,  one  or  more  neighbors  are  identified  for  each  particle,  and 
hence  evolution  equations  analogous  to  Eq.  (2)  can  be  constructed. 

This  model  has  been  implemented  and  tested  in  up  to  10  dimensions.  A  great  deal  of 
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effort  has  boon  concentrated  on  the  specification  of  (lie  coefficients  (analogous  to  in 

Rq.  2).  A  Gaussian  test  case  has  been  developed  which  sheds  light  on  the  specification. 

As  described  in  Section  3.  present  and  future  work  focuses  on  the  implementation  and 
testing  of  this  model. 


2.2  Stochastic  Lagrangian  Models 

The  mixing  model  described  in  the  previous  Section  fits  in  the  general  framework  of  stochastic 
Lagrangian  models.  Such  models  for  the  velocity,  composition  and  other  properties  following 
fluid  particles  can  be  used  to  effect  closure  of  the  Lagrangian  and  F.ulerian  pdf  equations. 
Hence  they  comprise  a  pdf  turbulence  model  applicable  to  turbulent  reactive  flows  such  as 
combustion. 

During  the  reporting  period,  two  substantial  papers  have  been  completed  on  the  topic  of 
stochastic  Lagrangian  models.  The  first.  Pope  ( 1 994),  is  a  tutorial  review  of  the  approach. 
The  second,  Pope  (1993),  explores  the  relationship  between  stochastic.  Lagrangian  models 
and  second-moment  closures. 

The  major  findings  and  contributions  of  the  latter  work  are  now  summarized.  To  every 
stochastic  Lagrangian  model  there  is  a  unique  corresponding  second-moment  closure.  In 
terms  of  the  second-order  tensor  that  defines  a  stochastic  Lagrangian  model,  correspond¬ 
ing  models  are  obtained  for  the  prcssurc-rate-of-strain  and  the  triple- velocity  correlations 
(that  appear  in  the  Reynolds-stress  equation)  and  for  the  pressure-scrambling  term  in  the 
scalar  flux  equation.  There  is  advantage  in  obtaining  second-moment  closures  via  this  route, 
because  the  resulting  models  automatically  guarantee  realizability. 

Some  new  stochastic  Lagrangian  models  are  presented  that  correspond  (either  exactly  or 
approximately)  to  popular  Revriolds-strcss  models. 
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2.3  Exact  Expression  for  Stationary  PDF 

Central  to  this  project,  and  indeed  to  many  issues  in  turbulence;,  is  the  shape  adopted  by 
the  pdf’s  of  different  flow  properties.  In  turbulent  combustion  the  importance  of  pdf’s  has 
been  recognized  for  over  twenty  years.  In  the  theoretical  turbulence  community  there  is  now 
much  interest  in  the  topic,  since  pdf’s  are  a  natural  diagnostic  for  intcrmittency. 

Together  with  Emily  Clung,  we  have  obtained  a  surprisingly  simple  exact  expression  for 
the  pdf  of  amplitude  of  any  stationary  random  process  (Pope  &  Ching,  1993).  Consider, 
for  example,  the  temperature  T(t )  measured  as  a  function  of  time  at  a  particular  point  in  a 
turbulent  reacting  flow.  Let  X (t)  be  the  standardization  of  T{t),  i.e. 


X(t)  =  (T(t)-(T))/nT: 


(5) 


where  (T)  and  aT  are  the  mean  and  standard  deviation  of  T{t).  The  result  we  have  obtained 
for  P{x),  the  pdf  of  X .  is 


P(x)  = 


C, 


(X*\x) 


exp 


{XW) 


(6) 


where  C\  is  a  normalization  constant,  and  (X|x)  denotes  the  expectation  of  the  second  time 
derivative  of  X(t)  conditional  on  X(t)  =  x. 

In  experiments,  simulations  and  modelling,  this  formula  provides  a  valuable  connection 
between  the  pdf  and  the  time  derivatives  of  the  signal.  As  diseased  by  Pope  &  Ching  (1993) 
it  sheds  light  on  the  tail-shape  of  pdf’s,  and  it  explains  the  success  of  a  previous  empirical 
expression  due  to  Ching. 


3  OUTLOOK 

The  principal  activity  in  the  final  year  of  the  project  will  be  to  test  and  implement  the  EMST 
model.  Preliminary  tests  (in  progress)  are  for  a  Gaussian  test  cause.  Then  the  model  will 
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be  tested  against  DNS  data  on  non-premixed  combustion,  obtained  at  Cornell  under  DOE 
funding. 

The  model  will  be  implemented  in  a  pdf  code,  combined  with  a  simplified  kinetics  scheme, 
and  applied  to  piloted  jet  diffusion  flames.  These  flames  provide  a  challenging  test  case.  Since 
blow-off  and  local  extinction  are  observed,  clearly  finite-rate  kinetics  are  important.  And  the 
application  of  existing  mixing  models  has  revealed  the  shortcomings  of  non-local  models. 


4  REFERENCES 

Chen,  H.,  Chen  S.,  and  Kraichnan,  R.H.  (1989)  Phys.  Rev.  Lett.,  63,  2657. 

Gao,  F.  (1991)  Phys.  Fluids  A,  3,  511. 

Kraichnan,  R.H.  (1990)  Phys.  Rev.  Lett.,  65,  575.  • 

Peters,  N.  (1984)  Prog.  Energy  Combust.  Sci.,  10,  319. 

Pope,  S.B.  (1982)  Combust.  Sci.  Technol.  28,  131. 

Pope,  S.B.  (1990)  Twenty-third  Symp.  (Int’l.)  on  Combust.,  The  Combustion  Institute,  p. 
591. 

Pope,  S.B.  (1991)  Theoretical  and  Computational  Fluid  Dynamics,  2,  255. 

Pope,  S.B.  (1993)  Phys.  Fluids  A  (submitted). 

Pope,  S.B.  (1994)  Annu.  Rev.  Fluid  Mech.  26  (to  be  published). 


5  PUBLICATIONS 

During  the  reporting  period,  the  following  papers  acknowledging  AFOSR  support  have  been 
submitted  for  publication. 

S.B.  Pope  (1994)  “Stochastic  Lagrangian  Models  for  Turbulence,”  Annual  Reviews  of  Fluid 


7 


Mechanics,  Vol.  26  (to  be  published). 

S.B.  Pope  (1993)  “On  the  Relationship  Between  Stochastic  Lagrangian  Models  of  Turbulence 
and  Second-Moment  Closures” 

S.B.  Pope  and  E.S.C.  Ching  (1992)  ‘Stationary  probability  density  functions  in  turbulence." 
Physics  of  Fluids  A  (to  be  published). 

S.  Subrarnanian  and  S.B.  Pope  (1992)  “Limitations  of  the  Amplitude  Mapping  Closure,” 
Bull.  Am.  Phys.  Soc.  37,  1805. 

6  PERSONNEL 

Professor  S.B.  Pope,  Principal  Investigator 
Tom  Dreeben,  Ph.D.  student 
Shankar  Subrarnanian  Ph.D.  student 
Dr.  Song  Fu,  Post-doc 

7  DEGREES  AWARDED 

None 

8  INTERACTIONS 

Talks  given  and  meetings  attended  during  the  reporting  period  are: 

•  Princeton  University,  Colloquium 

•  Stanford/NASA  Ames,  CTR  Advisory  Committee 

•  DOE  Contractors  Meeting,  Lake  Tahoe 


10%  AY,  1  month  summer 
50% 

50% 

50% 


8 


•  AFOSIl  Contractors  Meeting,  San  Diego 


•  21th  International  Combustion  Symposium.  Sydney,  Australia 

•  Stanford/NASA  Ames.  Center  for  Turbulence  Research:  Panel  Discussion  on  Turbulent 
Combustion;  Lectures  on  PDF  Methods 

•  General  Electric  CR&D,  Sandia  STAR  Meeting 

•  6th  Toyota  Conference,  .Japan 

•  Nagoya  University,  Japan.  Colloquium 

•  DOE  Review.  Bethesda.  Maryland 

•  APS/DFD  Meeting.  Tallahasse,  Florida 

•  Cornell  University,  Colloquium 

The  PI  is  involved  in  collaborative  work  with  researchers  at:  General  Electric.  Allison 
Gas  Turbine,  University  of  Stuttgart.  University  of  Sydney. 

9  INVENTIONS  AND  PATENTS 

None 


9 


